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ABSTRACT 

We study the formation of tidal tails in pairs of merging disk galaxies with structural 
properties motivated by current theories of cold dark matter (CDM) cosmologies. 
In a recent study, Dubinski, Mihos & Hernquist (1996) showed that the formation of 
prominent tidal tails can be strongly suppressed by massive and extended dark haloes. 
For the large halo-to-disk mass ratio expected in CDM cosmologies their sequence of 
modcls failed to produce strong tails like those observed in many well-known pairs of 
interacting galaxies. In order to test whether this effect can constrain the viability of 
CDM cosmologies, we construct N-body models of disk galaxies with structural prop- 
erties derived in analogy to the analytical work of Mo, Mao & White (1998). With a 
series of self-consistent collisionless simulations of galaxy-galaxy mergers we demon- 
strate that even the disks of very massive dark haloes have no problems developing 
long tidal tails, provided the halo spin parameter is large enough. We show that the 
halo-to-disk mass ratio is a poor indicator for the ability to produce tails. Instead, the 
relative size of disk and halo, or alternatively, the ratio of circular velocity to local 
escape speed at the half mass radius of the disk are more usehil criteria. This result 
holds in all CDM cosmologies. The length of tidal tails is thus unlikely to provide 
useful constraints on such models. 
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1 INTRODUCTION 

In standard hierarchical scenarios for galaxy formation, 
mergers of galaxies are common events that lead to the 
build-up of ever more massive galaxies. In fact, such merg- 
ers have been observed for a long time. There is now a large 
database of well studied examples of merging or strongly in- 
teracting disk galaxies, among the most prominent of them 
are NGC4038/39 (the Antennae), NGC4676 (the Mice), and 
NGC7252. Many of these pairs feature extended tidal tails, 
with a length that can reach more than 100 /i _1 kpc in pro- 
jection, or in the extreme case of IRAS 19254-7245 (the Su- 
perantennae) even ~ 305/i _1 kpc from tip to tip (Colina 
et al. 1991). 

The tails originate in close encounters of disk galax- 
ies, when the mutual tidal field ejects disk stars into arcing 
trajectories that lead to the formation of long tails point- 
ing way from the galaxies, and of bridges connecting them. 
This process was first demonstrated convincingly in a clas- 
sic paper by Toomre & Toomre (1972). Later White (1978, 
1979) computed the iirst fully selhconsistent 3-dimensional 
simulations of merging galaxies and established the rapidity 
of the orbital decay, and the structural resemblence of the 
merger remnants to elliptical galaxies. This work has been 
coniirmed and extended over the years by simulations with 
increasingly realistic initial conditions and ever better nu- 



merical resolution (Farouki & Shapiro 1982; Farouki et al. 
1983; Negroponte & White 1983; Barnes 1988, 1989, 1992; 
Hernquist 1992, 1993b; Barnes & Hernquist 1996). There 
have also been quite successful attempts to model particular 
interacting systems in detail, for example NGC7252 (Hib- 
bard & Mihos 1995) and NGC2442 (Mihos & Bothun 1997). 

Recently, Dubinski, Mihos & Hernquist (1996, hereafter 
DMH) studied the morphology of tidal tails in a series of 
merging models of disk galaxies with varying halo-to-disk 
mass ratio. In their sequence of four models, they kept the 
inner rotation curve very nearly constant and surrounded 
the disk and the bulge with ever more extended and mas- 
sive dark haloes. They found that with increasing mass of 
the dark halo, the resulting tidal tails became shorter and 
less massive. Their explanation for this effect is simple. For 
a fixed structure of the disk, a more massive halo leads to 
a deeper potential well and a higher encounter velocity. As 
a consequence, the duration and overall strength of the per- 
turbation to the disk is smaller, and the perturbed material 
cannot as easily climb out of the deeper potential well. 

Dubinski, Mihos, and Hernquist have followed up this 
study with an analysis of NGC7252 (Mihos et al. 1998), 
and an investigation of different dark matter proHles with 
a restricted 3-body code (Dubinski et al. 1997). Again they 
found that disk models with large halo-to-disk mass ratio 
were not able to produce prominent tidal tails. In particular, 
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they concluded that for mass ratios above 10:1 it should be 
exceedingly dirncult to make tails as long as those observed 
in systems like NGC7252 or NGC4038/39. Since the cur- 
rently favoured theoretical values of halo-to-disk mass are 
considerably higher than this, they speculated that there 
might be a connict with cold dark matter (CDM) cosmolo- 
gies. 

In this work we examine the tail-forming ability of re- 
alistic models of disk galaxies, where 'realistic' means that 
their structural properties are motivated to a large degree by 
current theories of CDM cosmologies. We derive the struc- 
tural properties of our disk galaxies according to the analytic 
model of Mo, Mao & White (1998, hereatter MMW), and 
we collide pairs of these galaxies in self-consistent N-body 
simulations. We adopt initial conditions for these merger 
simulations that are favourable for tail formation. 

We will demonstrate that the halo-to-disk mass ratio is 
not a particularly useful parameter for characterizing ability 
to make tidal tails. We find that it is the relative distribution 
of disk and halo material that is relevant, not the mass ratio 
itself. This conclusion was reached earlier by Barnes (1997, 
private communication) through analysis of a series of sim- 
ulations with halo/disk models differing both from those of 
DMH and from the CDM-based models we use here. 

We will show that realistic disk models in CDM cos- 
mologies have no problem producing long and massive tidal 
tails, provided the spin parameter of their dark halo is large 
enough. This result is practically independent of the adopted 
halo-to-disk mass ratio. 

This work is organized as follows. In Section 2 we de- 
scribe the structural properties of our disk models, and in 
Section 3 we discuss our techniques for setting up N-body 
representations of these models. A description of the dif- 
ferent simulations we have performed is given in Section 4, 
while Section 5 presents the results. Finally, we summarize 
and discuss our nndings in Section 6. 



2 MODELS OF DISK GALAXIES 

MMW have developed an analytical model for the structure 
of disk galaxies embedded in cold dark matter haloes. Their 
model rests on a number of simple yet plausible assump- 
tions, and it is very successful in reproducing the observed 
properties of disk galaxies. In particular, the predicted pop- 
ulation can match the slope and scatter of the Tully-Fisher 
relation as well as the properties of damped Lya absorbers 
in QSO spectra. We take their model as basis to derive the 
structural properties of our N-body models of disk galaxies. 
For denniteness, we briefly summarize the relevant assump- 
tions and equations. 



2.1 Dark haloes 

Using high-resolution N-body simulations, Navarro, Frenk & 
White (1996, 1997, hereatter NFW) established that haloes 
formed by the gravitational clustering of cold dark matter 
exhibit a universal structure. Suitably scaled, the density 
distribution of these dark matter haloes does not depend on 
cosmology. The NFW-profile is given by 

p(r) = pcrit. , . , ,f C —2 , (1) 



where p cr it. is the background density at the time of the halo 
formation, r s is a scale radius, and S c is a characteristic over- 
dcnsity. Note that the slope of this prohle is shallower than 
isothermal at the center, and it gradually steepens outward 
to an asymptotic slope of —3. Following NFW, we denne the 
virial radius r2oo as the radius with mean overdensity 200, 
i.e. it contains the virial mass 



4n 3 

M 20 o = 200p C rit. -g-r 20 o , 

and we dehne the concentration 

_ ^200 

r c 



(2) 



(3) 



of the halo. With these dennitions, the characteristic over- 
density is given by 



200 



3 ln(l + c) 
Further, let 



2 _ GM200 

«200 = 



1+c 



^200 



(4) 



(5) 



be the circular velocity at the virial radius. Given the con- 
centration c and the Hubble constant H(z), the radial den- 
sity pronle of a halo may then be specined by anyone of the 
parameters «200, »"200, or M2oo- In particular, we have 



M 2 



u 200 



WGH(z) 



and r 2 oo 



«200 
WH(z)' 



(6) 



2.2 Putting a disk into the halo 

We now put a stellar disk into an NFW halo according to 
the model of MMW. This rests on four key assumptions: 

(i) The mass M d of the disk is a given fraction m d of the 
halo mass. 

(ii) The spin J d of the disk is a given rraction j d of the 
angular momentum J of the halo. 

(iii) The disk has the structure of a thin exponential disk, 
and it is cold and centrifugally supported. 

(iv) Only disks that are dynamically stable against bar 
formation correspond to observable disk galaxies. 

The angular momentum J of a halo with total energy E is 
often characterized by the dimensionless spin parameter 



A : 



4__ 

GMi 



(7) 



According to N-body results (Warren et al. 1992; Lemson & 
Kauffmann 1998), the distribution of A is well approximated 
by 



p(A)dA = 



(2-kY^o 



exp 



(lnA-lnA) 2 



2a 2 



dA 
A 



(8) 



(r/r s ) (1 + r/r s ) 



with a = 0.5 and a typical value A = 0.05. This distribution 
is practically independent of cosmology, and of the mass and 
environment of the haloes (Lemson & Kauffmann 1998). The 
initial kinetic energy of the spherically symmetric halo may 
be computed by assuming that all particles move around the 
center on circular orbits, with speed equal to the circular 
velocity. This 'trick' results in 
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where 



GM 
2r 2 () 



/c, 



ln(l + c) 



[ln(l+c)-^]- 



(9) 



(10) 



Using the virial relation E = — -Ekin, the angular momentum 
of the halo then becomes 



(11) 



V /c ) 

We now put a fraction ma of the initial halo mass into 
a thin stellar disk with an exponential surface density, viz. 

R_ 

with E = M d /(2TvRl). Here M d = m d M 2 oo is the total 
mass of the disk and R d is its scale radius. The condition 



S(i?) = S exp(--^) 



(12) 



Jd = jd-J 



(13) 



will then determine the scale radius of the disk, because its 
spin is given by 

jd=Md r°°(^) 2 " (i?)Gxp (^) di? ' (M) 

where the circular velocity v c is the sum of two contribu- 
tions, namely 

V C (R) = R^ = Wc,disk(-R) + «c,dm(fl)- (15) 



2.3 Response of the dark matter pronle 

We take the gas, that later forms the disk, to be initially 
distributed just like the dark matter. However, the struc- 
ture of the dark halo will be changed when the disk forms in 
its center. We again follow MMW and assume that the dark 
matter reacts adiabatically to the disk formation. In par- 
ticular, we assume that the spherical symmetry of the halo 
is retained, and that the angular momentum of individual 
dark matter orbits is conserved. This latter condition may 
bc formulated as 



nM(n) = r{M f (r f ). 



(16) 



Here n and r f are the initial and final radii of some dark 
matter mass shell, M(r) gives the initial NFW mass profilc, 
and Mf (r) is the final cumulative mass profile after the disk 
is formed. M f (r) is the sum of the cumulative mass of the 
disk and the dark mass inside the initial radius, i.e. 



M f (r f ) = M d (r f ) + (1 - m d )M(n). 



(17) 



The tinal profile M h (r) of the dark matter halo is then given 
by 



M h (r) = M f (r) - M d (r). 



(18) 



For a given set of parameters «200, c, m d , j d , A, and 
a formation redshift z, the above equations uniquely deter- 
mine a disk model. Note that in practice the scale length 
R d of a disk needs to be determined iteratively in order to 
satisfy equations (13) and (17). 



2.4 Including a bulge 

In many galaxies, including the Milky Way, a central bulge 
population of stars is observed. For spirals like the Milky 
Way or of later type, the bulge mass is less than 20% of 
the disk mass. For this reason, the dynamical importance of 
the bulge in these systems should be small. However, there 
are also systems with a higher mass fraction in the bulge. 
While most of the models in this study do not have a bulge, 
we still want to investigate its possible inAuence on our re- 
sults. Hence we here generalize the above model to allow the 
option of a bulge. 

Bulges appear to be Aattened triaxial systems, that 
may be partly supported by rotation. However, Hernquist 
(1993b) found that it hardly matters for the density and 
velocity structure of merger remnants whether bulges are 
spinning or not. For simplicity, we therefore neglect a possi- 
ble Aattening of the bulges and model them as non-rotating 
spheroids with a spherical Hernquist profile of the form 



Pb(r) 



M h 



r h 



2tt r(r b + r) 3 ' 



(19) 



In analogy to the treatment of the disk, we assume that 
the bulge mass is a fraction m h of the halo mass. Since 
we take the bulge to be non-rotating it has lost its specific 
angular momentum either to the halo, or to the disk. We 
will assume that there is no angular momentum transport 
between the disk and the dark halo, and none between the 
disk and the bulge. In this case j d = m d . 

For simplicity, we further assume that the bulge scale 
radius r h is a fraction f h of that of the disk, i.e. r h = f h R d . 
Note that the disk half mass radius is 1.678 R d , while that 
of the bulge is 2.414 r h . 

It is then straight forward to generalize the above disk 
model to accomodate the bulge. The circular velocity of 
equation (15) gets an additional contribution from the bulge, 
i.e. 

Vc(R) = <disk(#) + v'l dm (R) + vl h (R), (20) 

with v cb (R) = GM h (R)/R. Further, equation (17) needs to 
be replaced by 

M f (r f ) = M d (r f ) + M h (r f ) + (1 - m d - m h )M(n), (21) 
and the dark mass protile of equation (18) now becomes 
M h (r) = M f (r) - M d (r) - M h (r). (22) 



3 N-BODY REALIZATIONS 

3.1 N-body realizations of model galaxies 

In order to construct near-equilibrium N-body realizations 
of our disk models, we need to initialize both positions and 
velocities of particles according to the solution of the col- 
lisionless Boltzmann equation (CBE)*. While the first can 
easily be done according to the derived mass distributions 
for halo, bulge, and disk, the latter is considerably more 
complicated. 

Instcad of attempting to solve the CBE directly, we 



also known as Vlasov cquation. 
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follow Hernquist (1993a) and assume that the velocity dis- 
tribution at a given point in space can be sumciently well 
approximated by a multivariate Gaussian. In this case, only 
the tirst two moments of the velocity distribution are needed. 
They can be obtained by taking moments of the CBE, a pro- 
cess that leads to a hierarchy of generalized Jeans equations 
(Magorrian & Binney 1994). 

For a static, axisymmetric system, the energy E and 
the angular momentum component L z are conserved along 
orbits. With the assumption that the distribution function 
depends only on E and L z one can show (Magorrian & Bin- 
ney 1994) that the first velocity moments are given by 



Vr = V z = VrV z = V z V^ = VrV^, = 0, 



-77 1 



<9$ 



Az p{R,z') — {R,z'), 



R d . -o-. 2 
^-p-dR ipV * )+V °> 



(23) 
(24) 

(25) 
(26) 

(27) 



where the azimuthal circular velocity is deiined as 

2 _ „9$ 
V *= R dR- 

Not specified by the Jeans equations is the azimuthal 
streaming v^, which can essentially be freely chosen in the 
context of the above approximations. This reAects the fact 
that the distribution hmction is even in L z ; the relative con- 
tribution of the parts with positive and negative L z can be 
arbitrarily chosen. 

We employ the assumption / = f{E,L z ) for the dark 
matter halo and the optional bulge. However, a realistic dis- 
tribution function for the disk has a more complicated struc- 
ture, and we will treat it slightly differently, as described 
below. 



3.1.1 Structure of the disk 

Real stellar disks have a iinite thickness. For their vertical 
structure we adopt the common choice of Spitzer's isother- 
mal sheet, viz. 

E( ^ech 2 (£). 



p d {R,z) 



2z 



(28) 



Here the thickness zo of the disk sets its 'temperature'. Most 
spiral galaxies seem to be consistent with a constant vertical 
scale length with a value of zo ~ 0.2Rd, which we will adopt 
in the following. 

The distribution function of the disk depends on more 
than just two conserved quantities, hence it is unrealistic 
to assume an isotropic velocity dispersion. However, we will 
keep the assumption that the velocity ellipsoid is aligned 
with the coordinate axes. Then equation (25) remains valid, 
and we use it to compute a z {R,z). Note that due to the 
radial variation of E(i?) and the presence of the halo, the 
vertical velocity structure of the disk will not be exactly 
isothermal. 

We further employ the epicycle approximation (Binney 
& Tremaine 1987, chapter 4) to relate the radial and az- 
imuthal velocity dispersions by 



„2 

2 °r 



(29) 



Here we have deiined 
4 9$ 



7 



k?R dR' 



(30) 



(31) 



and the epicyclic frequency k as 

2 = _. __ _!_. 
K = R dR + dR 2 ' 

The epicycle approximation also implies that the asymmet- 
ric drift is small. We neglect it altogether and set the stream- 
ing velocity to be equal to the circular velocity, i.e. v$ = v c . 
For simplicity, we also continue to assume that o R = o~\. 
With these assumptions the velocity structure of the disk is 
fully specified. 

Note that the thickness of the disk must be chosen large 
enough to fulfil Toomre's stability criterion, which requires 



Q 



o-rk 



> 1 



(32) 



3.36GS 

to ensure local stability in differentially rotating disks. 

For our models, the minimum value of Q is about 1.4. 
Hence, we could have taken somewhat colder disks, which 
might produce sharper tidal tails. 

3.1.2 Rotation of the halo 

For consistency, we want to properly represent the angular 
momentum carried by the dark matter, although we do not 
expect that halo rotation will have a strong innuence on tail 
formation. 

We model the streaming velocity of the dark halo as 
some fixed haction / s of the local azimuthal circular velocity, 
i.e. 



' = fsVc 



(33) 



If the specific angular momentum of the dark matter is con- 
served during disk formation, the factor / s stays fixed as 
well. Hence it can be computed for the initial NFW-profile. 
For the streaming of equation (33) the initial angular mo- 
mentum is 



J 



[Ml + c)-^] 
where g c is the integral 



Jo 



ln(l + x) - 



l+x 



3 

12 



(1 



-d:r. 



(34) 



(35) 



Comparing this with equation (11) we see that / s is given 
by 



ln(l + c) 



1 + c 



(36) 



The quantity / S /A varies only weakly, e.g. it takes the value 
4.5 for c = 1, and reduces to 3.3 for c = 100. 



3.1.3 Halo truncation 

For the halo we encounter a slight technical problem, since 
the cumulative mass distribution of the NFW profile actu- 
ally diverges for large r. This is simply due to the fact that 
the NFW-profile in the form of equation (1) is not valid out 
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to arbitrarily large distances; it just provides a good fit to 
the pronle up to about the virial radius. 

Instead of truncating the prorile sharply at the virial 
radius we rather want to derive N-body models where the 
density fades out smoothly. For simplicity, we have chosen 
an exponential cut-off that sets in at the virial radius and 
turns off the prorile on a scale r s , viz. 



p(r) = 



Pcrit.^c 

z(\ + c) 2 



'•21)11 



exp 



r - rgoo \ 



(37) 



for r > r 2 oo- The power law exponent a allows a smooth 
transition of the prorile at r 20 o- We select a such that the 
logarithmic slope 



d . 
r — lnp 

dr 



1 + 3c 
' \ + c 



(38) 



of the prorile at the virial radius is continuous. This implies 
a = c + n. 

Note that this truncation results in some additional 
halo mass beyond the virial radius. The total mass M to t 
is roughly 10% larger than M 2 oo- However, we want to keep 
our dehnitions of disk and bulge masses in terms of the virial 
mass. As a consequence we need to slightly modiry equation 
(17). It becomes 

M 200 , 



M f (r f ) = M d (r { ) + M b (r f ) + 1 - 



M u 



-(m d + m b ) 



M(n 



3.1.4 Numerical procedure 

Finally, we briefly describe our computer code to set up a 
galaxy according to the above model. The following steps 
are followed: 

(i) Particle positions are initialized according to the den- 
sity prohles for halo, disk, and bulge. 

(ii) We then compute the velocity dispersions on a nne 
logarithmic mesh in the (R, z)-plane. For this purpose, we 
compute the integrals 



dz' p 



9$ 

6V 



(39) 



numerically at the grid points. This also determines v\ — v1. 
Note that the density refers only to the component under 
consideration, while the potential is given by all the mate- 
rial. For the halo and the bulge, we use equation (26) and 
find v^ by numerically differentiating pvj^ in this plane. For 

the disk, we use the epicycle approximation to determine cr^, 
and we set the azimuthal streaming equal to the local circu- 
lar velocity. The streaming of the halo is given by equation 
(33), while that of the bulge (if present) is set to zero. 

It should be noted that the integrals of equation (39) re- 
quire elaborate numerical techniques, since the computation 
of the combined force field is nontrivial. 

(iii) Finally, particle velocities are initialized by drawing 
random numbers from multivariate Gaussians with disper- 
sions interpolated from the (R, z)-grid to the particle posi- 
tions. 

This scheme is similar to that of Hernquist (1993a), 
although our numerical procedure and treatment of the disk 
is somewhat different. 



Table 1. Paramctcrs of our basic set of six disk models. All 
of them have the same total mass corresponding to V 200 = 
160 kms -1 , and an initial halo concentration of c = 15. Prom 
the possible combinations of X,m d G {0.025,0.5,0.1} wc con- 
sider only those modcls that are stable according to the criterion 
A > m d . We assume that the disk material conserves its specinc 
angular momcntum, i.c. j d = m d . The two tables illustrate our 
labeling of the models and the resulting disk scale lengths. Note 
that these models have no bulgc. 
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Figure 2. Sizcs of disks. The two solid lincs show the scale length 
R d of disks as a mnction of the spin parameter A. The solid lincs 
are for haloes with c = 15, whilc thc dottcd curves refer to c = 5. 
In both cases, the heavy curves are for m d = 0.025, and the thin 
ones for m d = 0.1. The dashed line shows the distribution p(A) 
of A cxpcctcd in CDM cosmologies. 



4 SIMULATIONS 
4.1 Models 

4-1.1 The basic disk models 

We have constructed a basic set of six disk models with 
a constant total mass corresponding to V 20 o = 160 km s _1 
and a concentration of c = 15. These models are those com- 
binations of A, m d £ {0.025,0.5,0.1} that result in stable 
cold disks, i.e. that have A > m d (the stability criterion is 
discussed in more detail below). We label these models A 
to F, as outlined in Table 1. Also given in this Table are 
the resulting disk scale lengths. Note that we always assume 
jd = rn d , i.e. the specihc angular momentum content of the 
disk material is exactly conserved during disk formation. In 
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m d = 0.100 X = 0.100 



2 3 

R/R d 

Figure 1. Rotation curves of disk modcls A to F. The top three panels show the full rotation curve of models A-C out to the virial 
radius. In the other six panels, the radial coordinatc is normalized to the disk scale length, and we plot the inner rotation curve out to 
5 disk scale lengths. In each panel, the dotted curve gives the contribution of the dark matter, the dashcd that of the disk, and the solid 
line is the total rotation curve. In all the models, the total mass of the galaxies is equal, and corresponds to an initial NFW dark matter 
pronle with «200 = 160 km s _1 and c = 15. 



contrast to this, gas-dynamical simulations of disk forma- 
tion (Navarro & White 1994; Navarro & Steinmetz 1997; 
Weil et al. 1998) have typically led to a loss of angular mo- 
mentum from the gas to the halo. As a result, the disks 
formed in these simulations were much too small to be iden- 
tihed with real spiral galaxies. However, as all of the above 
authors note, this angular momentum problem may well be 
due to an insufncient treatment of feedback processes. 



Note that in the model of MMW, the structure of the 
disk galaxies depends only on A' = (j d /m d )A. Hence, angu- 
lar momentum loss from the disk (j d < m d ) has the same 
effect as lowering the value of A. 



4-1.2 Rotation curues 

In Figure 1 we show the rotation curves of our six primary 
models described in Table 1. For each model, we give the 
inner rotation curve out to 5 disk scale lengths, which is 
about the accessible regime in most disk galaxies. For models 
A to C, we also show the full rotation curve out to the virial 
radius of 160 fe _1 kpc. 

Several interesting trends may be observed. In the mod- 
els A, B, and C, only the spin parameter is increased. This 
leads to larger disks with roughly i? d tx A. The dependence 
of Rd on the spin parameter A is shown in Figure 2. However, 
the smaller disks pull in the dark matter more strongly, lead- 
ing to a larger concentration of the dark matter for smaller 
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Table 2. List of runs. The table gives the orbital angular mo- 
mcntum of thc different runs in terms of the minimum Kcplcrian 
separation i?kep- All the runs had an initial galaxy separation of 
-Rstart = 320h. _1 kpc, and were set-up on a parabolic encounter 
with zero total energy. Each of the runs is a collision between 
idcntical disk modcls. The latter is specihed by the initial char- 
acter of the labels. The bottom table gives the actual separation 
Rmin of thc disks in thcir nrst cncounter. 







^kep 


/i-^kpc] 




3.5 


7.0 


14.0 


28.0 56.0 


112.0 


AO 


Al 


A2 






BO 


Bl 


B2 






CO 


Cl 


C2 


C3 C4 


C5 






C2r 










C2i 








Dl 












E2 










F2 







Tl 
Ul 
VI 
Wl 







-^min 


/i-^kpc] 






A0 


6.4 


B2 


23.7 


Dl 


10.8 


B0 


8.3 


C2 


22.4 


E2 


18.2 


CO 


7.4 


C2r 


19.4 


F2 


22.4 


Al 


10.5 


C2i 


21.2 


Wl 


10.9 


Bl 


14.0 


C3 


37.4 


Tl 


8.3 


Cl 


13.6 


C4 


58.2 


Ul 


9.4 


A2 


21.3 


C5 


109.6 


VI 


10.8 



disks. This effect reduces the differences between the rota- 
tion curves when their radial coordinate is normalized to the 
disk scale length. 

On the other hand, for very small disk mass, the dark 
matter prohle will be nearly unaffected by the disk forma- 
tion. In this limit, the disk stars behave more or less like test 
particles in the dark matter potential, yet the size of the disk 
is still determined by the halo spin parameter. Also, in this 
limiting case of a massless disk it is quite clear that the 
mass ratio between halo and disk must be irrelevant for the 
formation of stellar tidal tails. 

4-1.3 A model with a masswe bulge 

We also consider a model with a massive bulge, designed to 
have a similar rotation curve to DMH's models, and hence 
being more directly comparable to them than our standard 
models. We also adopt their relatively high disk-to-bulge 
mass ratio of 2:1. In detail, our parameters for this model, 
which we call 'W', are c = 15, W200 = 160 km s _1 , A = 0.05, 
m d = 2/3 x 0.05, m b = 1/3 x 0.05, and / b = 0.1. This 
results in the rotation curves shown in Figure 3. While the 
inner rotation curve is very similar to that of DMH, the 
contribution of the dark matter is much more important in 
our model, even at small radii. In particular, it is always 
much larger than the contribution of the disk. 
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Figure 3. Rotation curves of the bulge model W. The top panel 
shows the full rotation curve, and the lower panel displays the 
inner rotation curve out to 5 disk scale lengths. The dashed line 
is the contribution of the disk, thc dot-dashcd that of the bulgc, 
the dotted that of the dark halo, and the solid line gives the total 
curve. 

4-1-4 The amount of dark mass in the disks 

As the rotation curves in Fig. 1 show, all our 'basic' disk 
models A-F are gravitationally dominated by dark matter, 
even in the innermost regions of the disks. While the pres- 
ence of a dark matter halo has been comdncingly demon- 
strated by the natness of observed rotation curves, there is 
still a controversy about the amount of dark matter in the 
inner regions of disk galaxies. 

This controversy has arisen, because the decomposition 
of an observed rotation curve into a stellar and a dark mat- 
ter component is rather ambiguous, since the result depends 
strongly on assumptions about the dark matter profile and 
the mass-to-light ratio (Navarro 1998). Traditionally, rota- 
tion curves have therefore been fitted using the 'maximum- 
disk' hypothesis, i.e. the largest mass possible is assigned to 
the disk consistent with the rotation curve. The recent work 
of Debattista & Sellwood (1998) suggests that the central 
density of dark matter in barred galaxies should be low, thus 
supporting the maximum-disk hypothesis. However, others 
(e.g. van der Kruit 1995) maintain that the contribution 
of dark matter in the inner regions of disk galaxies must 
be substantial. There is also recent observational work that 
supports this assertion (Quillen & Sarajedini 1997). We also 
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Figure 4. Contours of the disk stability paramctcr e m (solid 
lines) for c = 15. According to Syer, Mao & Mo Syer et al. (1998) 
disks with e m > 0.75 should be stable against bar formation, whilc 
the earlier work of Efstathiou, Lakc & Ncgropontc Efstathiou 
et al. (1982) gives the condition e m > 1.1. Disk modcls lying 
below the thick dashed line are dominated by the disk gravity 
at the maximum of the disk rotation curve, thus the hatched 
rcgion shows thc parameter space that may contain stable disks 
which are not everywhere dominated by dark matter. Two of our 
modcls, 'T' and 'U' (stars), lie in this region. The other models 
are indicated as triangles. 



note that even if the inner rotation curve of a galaxy can be 
well accounted for by a disk component alone, this does not 
provide evidence for the absence of dark matter in the inner 
parts of the disk. 

The theoretical results employed in this work predict 
dark matter pronles with a central density cusp, and a re- 
sulting strong contribution of dark matter in the inner disks. 
We now brieny discuss, to what extent the model of MMW 
may also accommodate galaxies that are maximum-disk, or 
at least somewhat closer to it. In order to make the self- 
gravity of the disk more important, we can either reduce 
A (making the disk smaller), increase m d (making the disk 
heavier), or lower c (reducing the concentration of the halo). 
However, for fixed m d and fixed c, disk-stability poses a 
lower limit on A. Thin, fully self-gravitating disks have been 
known to be violently bar-unstable for a long time, a fact 
that suggested to Ostriker & Peebles (1973) that there must 
be dark matter that stabilizes the disks. Later, Efstathiou, 
Lake & Negroponte (1982) used N-body simulations to de- 
rive the stability criterion e m > 1.1 for the disk, where 



400 



(40) 



and V max 1S the maximum rotation velocity. Recently, Syer, 
Mao & Mo (1998) connrmed that e m is a good diagnostic 
for bar-instability, although they found a somewhat weaker 
stability criterion, e m > 0.75. 

In Figure 4 we show contours of e m in the m d -A plane. 
Also shown is the region, where the disk gravity at the max- 
imum of the disk rotation curve is larger than the contri- 
bution by the dark matter. There is thus a small region 
of parameter space (hatched) where the galaxies are disk- 
dominated, but where they should be still stable against 
bar formation according to Syer et al. (1998). Incidentally, 
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Figure 5. Inner rotation curves of thc modcls T, U, and V out 
to 5 disk scale lengths. The dashcd line is the contribution of the 
disk, the dotted that of the dark halo, and the solid line gives the 
total rotation curve. 



e m > 1.0 corresponds very closely to the condition A > m d , 
the choice we employed so far. 

Note that lowering A may not only render a disk un- 
stable, it will also make it substantially smaller. Hence, the 
observed sizes of disk galaxies can also provide a lower bound 
on A. Indeed, the results of MMW suggest that disks become 
too small if they lose a substantial part of their angular mo- 
mentum to the dark halo. Since these arguments disfavour 
low A, one may rather try to increase m d to make the grav- 
ity of the disk more important. In principle, we expect that 
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the universal cosmic baryon fraction poses an upper limit 
on m<j, while the actual value of rrid could be a lot smaller 
if the efnciency of disk formation is low. Taking a big bang 
nucleosynthesis value of J7b — 0.015/j -2 (Copi et al. 1995) 
for the baryon density, rrid should be smaller than 0.06 in a 
critical density universe with a Hubble constant of h — 0.5. 
However, clusters of galaxies suggest that the baryon frac- 
tion is larger by at least a factor of three (White et al. 1993). 
Note that in a low density universe this can be reconciled 
with cosmic nucleosynthesis. If flo is as low as 0.2, the limit 
on m<j goes up to about 0.15-0.2. 

We examine these possibilities to a limited extent with 
three additional models which are disk-dominated in the in- 
ner regions. It is interesting to see how their tidal tails fit into 
the systematic properties of the other models. We label one 
of these models 'T', and give it the parameters rrid = 0.1, 
A = 0.05, c = 15, «2oo = 160kms _1 , and ra D = 0. For 
a further model, called 'U', we instead adopt rrid = 0.15, 
A = 0.08, i.e. here we make the disk substantially more mas- 
sive. Finally, we consider a model 'V with a smaller con- 
centration of the halo. Here we use c = 5, m d = 0.1, and 
A = 0.05. Hence this model is consistent with the value of 
c ~ 5 favoured by Navarro (1998) in a recent analysis of a 
sample of spiral galaxies. Note that such low concentrations 
are theoretically expected for flat, low-density universes. 

The inner rotation curves of these models are shown in 
Fig. 5. The stability parameter for them is e m ~ 0.84. Hence 
they lie in the hatched region of Figure 4. In contrast to 
the other models, we here chose o\r = 2.0 a z for the velocity 
structure of the disk to prevent it developing a bar before the 
galaxies collide. This raises Q to about 2.0, while it would 
have been Q ~ 1.0 for our conventional choice for or. 

4.2 Collision simulations 

The most favourable condition for making tidal tails are 
prograde encounters where the spin vectors of the disks are 
aligned with the orbital angular momentum. In this situa- 
tion, the approximate resonance between the disk rotation 
and the orbital angular frequency amplines the perturbation 
of particlc orbits on the far sides of the disks, since they stay 
for a longer time in the region of the strongest tidal field. 

Since we here try to achieve as prominent tails as pos- 
sible, we usually set up our disk-disk collisions on prograde 
parabolic orbits. For simplicity we run only symmetric en- 
counters between pairs of identical models; that is we collide 
model A with A, B with B, and so on. 

We always chose the initial separation of the galaxies 
to be twice the virial radius. i.e. i? s tart = 320ft -1 kpc. The 
remaining undetermined orbital parameter is the orbital an- 
gular momentum. We specity it in terms of the minimum 
separation 7?k cp the galaxies would reach if they were point 
masses moving on the corresponding Keplerian orbit. In re- 
ality, once the galaxies overlap they will start to deviate 
from this trajectory due to dynamical friction. As a conse- 
quence, the actual separation i? m i n of the galaxies in their 
first encounter will generally be larger than i?k cp - 

We examined three main choices for i?k cp , 3.5, 7 and 
14/i -1 kpc. For each of the models A, B, and C, we have 
run all three of these combinations, while we restricted our- 
selves to just one 'impact parameter' for the other models. 
Additionally, we simulated a set of wider encounters with 



Rkcp = 28, 56, 112/i -1 kpc for the C model. Runs labeled 
'A0', 'B0', etc. refer to Rk cp = 3.5/i -1 kpc, those containing 
the digits 1 or 2 to 7/i -1 kpc and 14 h~ kpc, respectively. 

We also simulated two additional versions of run C2 
where the disks do not have a prograde orientation. In the 
collision 'C2r' both disks are retrograde, i.e. their spins are 
just Aipped, while in the model 'C2i' they are inclined by 
90° relative to the orbital plane. 

With respect to the bulge model and the disk- 
dominated models, we have only run one simulation in 
each case (Wl', 'Tl', 'Ul', 'VI'), an encounter at i?k cp = 
7/i -1 kpc. Table 2 gives an overview of all these runs. Also 
shown in this table are the actual minimum separations i? m i n 
of the centers of the disks in their first encounter. We here 
denned the center of a disk as its densest point, and used 
a kernel interpolation like in smoothed particle hydrody- 
namics (SPH) to estimate the density of particles. However, 
using simply the center-of-mass of the indiviual disks gives 
similar results. 

When the galaxies start to overlap, the interaction po- 
tential between them becomes shallower than that of the 
corresponding point masses. This effect will make the orbits 
wider than the Keplerian expectation. However, the galaxies 
are also slowed down by dynamical friction, an effect that 
brings the galaxies closer together. Both mechanisms com- 
pete with each other, and their relative strength depends on 
the distribution of mass inside the galaxy. As Table 2 shows, 
the minimum separations i? m i n are usually somewhat larger 
than the corresponding Keplerian value i?k cp . Note however, 
that the measurement of i? m i n has an uncertainty of order 
±1 /i -1 kpc, because we stored only a limited number of out- 
put times. 

4.3 Numerical techniques 

All the simulations in this work have been run with 
our GADGET-code (GAlaxies with Dark matter and Gas 
intEracT). It is a newly written SPH-code in C, specincally 
designed for the simulation of galaxy formation and interac- 
tion problems. The gravitational interaction is either com- 
puted with the special-purpose hardware GRAPE (if avail- 
able) or with a TREE code. Time integration is performed 
with a multi-timelevel leapfrog intcgrator, with the ther- 
mal energy equation being integrated semi-implicitly. Time- 
critical routines (i.e. the force computation and neighbour 
search with the TREE or GRAPE) have been pronled and 
optimized extensively. 

In this work, the SPH part of GADGET is not used; we 
just treat dark matter and stellar material as collisionless 
particles. Futher details of GADGET and an application of 
its SPH-capabilities will be described in future work. 

For all of the basic models A to F, and for T to V, 
we used 20000 particles to represent each disk, and 30000 
particles for each halo, hence each simulation had a total of 
100000 particles. We chose a gravitational sottening length of 
0.4 /i -1 kpc for the dark matter, and 0.1 /i -1 kpc for the disk. 
Time integration was performed with high enough accuracy, 
such that the total energy was conserved to better than 0.8% 
in all runs. 

For the bulge model W, we used an additional 10000 
particles for each bulge, and we employed a softening of 
0.1 /i -1 kpc for the bulge as well. 
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Some of the models have been integrated using GRAPE 
(A0, B0, B2, E2, F2, Tl, Ul, Wl), the others with the TREE 
code. For the latter we used the cell opening criterion of 
Dubinski (1996) with 9 = 1.0, we included quadrupole mo- 
ments, and we matched the spline softening of the TREE 
code to the Plummer softening of GRAPE cited above. 

Each simulation was run for 2.6 internal time units, or 
0.26 Hubble times, corresponding to 2.54 x 10 9 /i _1 yr. At 
this point of time, the merger remnants are not yet fully 
relaxed, but the tidal tails have already largely decayed. 



5 RESULTS 

5.1 Dynamical evolution of the models 

Figure 6 shows a representative example of the time evo- 
lution of one of our runs (Bl). Overall, the models follow 
the well-known behaviour of close encounters of pairs of 
disk galaxies. When the galaxies reach orbital centre, vio- 
lent tidal forces induce a bar instability that quickly trans- 
forms the disks into a pair of open bisymmetric spirals. 
Simultancously, disk material from the far side of the en- 
counter is ejected by the tidal field into arcing trajectories 
that later form tidal arms. Material from the near side is 
drawn towards the companion, giving rise to bridges between 
the galaxies as they temporarily separate again. While the 
bridges are destroyed when the galaxies come back together 
for a second time, the tails can survive and grow for a longer 
time in the relatively quiet regions of the outer potential. 

Nevertheless, the dynamical evolution of the tidal tails 
is quite rapid. After their initial phase of expansion, the 
most strongly bound material in the inner region of the tail 
quickly starts to rain back onto the merging pair. Eventually, 
this also happens to material progressively further out, such 
that the surface density and prominence of the tidal tails 
quickly decrease with time. 



5.2 Comparison of tidal tails 

Depending on the strength of the tidal response, the tidal 
tails can contain a varying amount of mass, and reach differ- 
ent lengths. By comparing the time evolution of the models 
A to C, such differences are readily apparent. For example, 
the tails of the larger disks of run Cl are much more mas- 
sive and prominent than those of run Bl, while the tails of 
the small disks of simulation Al are rather thin and anemic. 
However, the spatial extent of the tails is quite compara- 
ble in the models. When normalized to the initial disk scale 
length the anemic tails of the A-models are even longer than 
those of the C-models. 

These trends are clearly visible in Figures 7 and 8, where 
we compare different runs at the same time, approximately 
corresponding to the moment when the tails are most im- 
pressive. Note however, that due to the rapid evolution of 
the morphology of the tidal tails it is not easy to comparc 
different models at exactly equivalent times of dynamical 
cvolution. 

When models with different impact parameters are 
compared, some finer trends in the tail morphology may be 
observed. With growing impact parameter (labels — » 2), 
the bridges between the galaxies become more pronounced, 



and the tails are slightly more curved. This is related to the 
larger orbital angular momentum of these encounters. 

Note that the runs A, B, and C of Figure 7 all have a 
large halo-to-disk mass ratio of more than 40:1, i.e. according 
to DMH none of these models should have produced promi- 
nent tails. Nevertheless, the spatial extent of the tails is of- 
ten quite large. For example, the B-models with a disk scale 
length of Rd ~ 3.5 ft _1 kpc produce tails reaching 200 /i _1 kpc 
in length, i.e. about 60 times the original disk scale length. 

When tails of models with different halo-to-disk mass 
ratio but equal spin parameter are compared, it becomes 
clear that the mass ratio cannot be the relevant parameter 
that decides whether tails form or not. For example, the 
tails of runs E2 and F2 in Figure 7 may be compared to the 
ones of run C2. Despite a variation of the mass ratio by a 
factor 4 or so, the tails are almost equally strong in these 
three simulations. This is clearly due to the approximately 
equal size of the disks in these models. Because the dark 
matter is gravitationally dominant even in the regions of 
the disks, the disk stars behave almost like test particles in 
the gravitational potential of the dark halo. In this limiting 
case it is clear, that only the location of the disk material 
inside the dark halo determines the disk response, i.e. it is 
the relative size of disk and dark halo that matters. 



5.3 An indicator for tidal response 

As we have seen above, knowledge of the halo-to-disk mass 
ratio is by no means sumcient to predict how prone a par- 
ticular galaxy model is to tail formation. MMW suggested 
using the quantity 



£ = 



v e (R) 



v c (R) 



(41) 



as a more suitable indicator. £ compares the depth of the 
potential well with the specinc kinetic energy of the disk ma- 
terial. The quantity £ also arises, when one tries to estimate 
the relative increase AE of specinc kinetic energy of disk 
stars in a nearly head-on encounter of identical disk models. 
In this situation one nnds (Binney & Tremaine 1987; Mo, 
Mao & White 1998) 



AE/v 2 c 



(42) 



This suggests the use of £ as an indicator for the susceptibil- 
ity of a disk model to tidal perturbations. In order to obtain 
a typical value for £ we evaluate it at 7? = 27?d, which is 
about the half mass radius of the disk. 

We now want to test how well £ works as an indicator 
for the ability of a particular disk-halo model to develop 
massive tails. Here one encounters two immediate problems. 

First, the tidal response of a disk depends on the or- 
bital parameters of the encounter with its companion. For 
cxample, if a disk is tilted against the orbital plane, the tidal 
forces felt by the bulk of the disk material will generally be 
smaller, resulting in a less prominent tidal tail. Similarly, a 
change of the impact parameter or the orbital energy can af- 
fect the tidal response. We here do not intend to investigate 
the complete parameter space. Rather, we focus on colli- 
sions that produce the strongest tails possible for mergers 
of a given disk model. 
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Figure 6. Time evolution of run Bl. The panels show the disk particles projected onto the orbital plane. The length units labcling thc 
axes are given in h~ kpc. The elapsed time since the start of the simulation (upper left corners) is given in units of 0.1 Hubblc times, 
or 9.8 x 10 s h _1 yr. 
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Figure 7. Comparison of tidal tails between different runs. All the panels show the disk particles projected onto the orbital plane at 
the same time since the start of the simulation. The top 9 panels display runs that involve the disk modcls A-C for the set of impact 
parameters R^ ep = 3.5, 7, and 14 /i — x kpc (indicatcd by the digits 0, 1, and 2, respectively, in the labels of the runs). From A to C, the 
spin parameter A increases in the sequence 0.025, 0.05, and 0.1, but all three models have an equal disk mass given by = 0.025. 
In the lower three panels we show additional collisons of model C with wider impact paramctcrs in thc sequence R^ ep = 28, 56, and 
112 Ji-^kpc. 
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Figure 8. Comparison of tidal tails between different runs. The panels show the disk particlcs projccted onto the orbital plane at the 
same time since the start of the simulation. Runs Dl and E2 involve models with twice as hcavy disks (m d = 0.05) than A-C, whilc 
F2 has m d = 0.1. The D-model has a spin paramctcr of A = 0.05, whilc the E- and F-disks have A = 0.1. The T (m d = 0.1, A = 0.05) 
and U (m d = 0.15, A = 0.08) modcls arc disk-dominatcd in their inner regions, while model W (m d = 0.034, m^ = 0.016 A = 0.05) also 
contains a bulgc. Thc digit in the label of the runs parameterizes the impact parameter; 0, 1, and 2 are for R^ cp = 3.5, 7, or 14/i _1 kpc, 
rcspcctively. 
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Figure 9. Tidal rcsponse T as a function of time for three selected 
runs. We define T as thc mass rraction of cach disk that reaches a 
distance of more than 10R d to its center-of-mass. Also indicatcd 
as horizontal lines are the values of T e g that we take as measure 
for the tidal response of the disk. Since we only stored simulation 
outputs with a spacing of 0.1 time units, the measurement of T e ff 
is slightly uncertain. 



In this spirit, prograde encounters are an obvious choice 
for the orientations of the disks, since this conhguration 
has repeatedly been shown to produce the strongest tails, 
and we here conhrm this with the runs C2r and C2i. We 
use parabolic orbits, because they are plausible candidates 
for real interacting galaxies if one assumes that they are 
coming together for their first time. However, even when 
this assumption is dropped, DMH showed that moderately 
bound orbits are no more effective in producing tails than 
zero energy orbits. With respect to the 'impact' parame- 
ter 7?kep, we have examined a range of different choices and 
found that the tails appear to be of maximum strength for 
Rmin ~ 2 — 6Rd, i.e. in collisions where the disks pass each 
other at a distance of a few disk scale lengths. Hence, our 
merger simulations of equal disk galaxies have been set up 
to exhibit the most favourable conditions for tail formation, 
and should indeed produce the strongest tails possible for 
these disk models. 

Second, it is not obvious how to define the mass or 
extent of the tidal tails in an objective way. This is hirther 
complicated by the rapid dynamical evolution of the tails, 
which makes it difficult to compare simulations that may 
form their tails at different times. 

In ordcr to solve this problem and measure the strength 
of the tidal response, we have come up with the following 
scheme. We start by defining the quantity T to be the mass 
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Figure 10. Tidal response T e g of thc diffcrcnt runs vcrsus thc value of £ for the corresponding disk model, where £ is cvaluated at 
R = 2i?d- The measurement of T e ff is slightly uncertain, and T e ff also dcpcnds on dctails of thc orbital paramcters of the encounter. 
However, our runs are designed to produce very nearly the strongest tails possible for a given disk modcl. At fixcd £, it should thcrefore 
be hard to nnd a modcl that gives a higher value for T e ff than the maximum of our runs. In other words, there should be no modcl in 
the region above the indicated dashed line. 
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Figure 11. Tidal response T e ff of the different runs versus the disk-to-halo mass ratio (left panel), and versus M(2R^)/M (right pancl). 
The quantity M(2R^)/M is a measure of the fraction of the total mass in the disk rcgion. 



fraction of each disk that reaches a distance of more than 
10-Rd from its center-of-mass, where Rd is the original scale 
length of the unperturbed disk. 

In Figure 9 we show examples for the time evolution of 
T. Shortly after the disks come together for the first time, 
T jumps up, reaches a maximum, and slowly decays, until 
the disks are scrambled up in their second encounter and T 
loses its initial meaning. Note that the different runs reach 
their maximum of T at different times. In order to compare 



them on an equal footing, we therefore define an effective 
response T e g as the peak value reached by T. 

In Figure 10 we plot the tidal response T e ff of our runs 
versus the value of £, evaluated at R = 2R d . Although our 
coverage of parameter space is limited, it is nevertheless clear 
that there is a correlation between T e ff and £. Of course, 
there is some uncertainty in the measurement of T e ff and this 
introduces some scatter. Also, for fixed £, the tidal response 
T e ff depends somewhat on the impact parameter _Rkc P , and 
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Figure 12. Tidal tails of simulations C2r and C2i. The panels show the disk particles, projected either onto the xy-plane (orbital plane), 
or the yz-plane. Apart from a different orientation of the disks, runs C2r and C2i are identical to C2. In C2r, both disks arc rctrogradc, 
while in C2i the spins of the disks lie in the orbital planc, pointing along thc x-axis. These runs may be compared with C2 shown in 
Figure 7. 



it might also have a slight dependence on m d . However, to 
the extent that our simulations really produce the strongest 
tails possible for our disk models, Fig. 10 shows that it will 
be exceedingly hard to find a model that produces tails that 
lie above the dashed diagonal line. This establishes that £ is 
a good indicator for the maximum tidal response T e g, that 
may be obtained for a given class of disk models. In par- 
ticular, models with £ > 8 should be unable to produce 
strong tails. This is in excellent agreement with the analysis 
of MMW, who estimated £ = 4.2, 5.5, 7.2, and 9.3 (in order 
of increasing halo mass) for the sequence of four models of 
DMH; in agreement with Figure 10, the last two of these 
models failed to produce prominent tails. 

Figure 10 may also be compared to the two panels of 
Figure 11. In the left panel, we plot the tidal response ver- 
sus the disk-to-halo mass ratio. This again shows, that the 
disk-to-halo mass ratio is not a good indicator for the abil- 
ity to form tidal tails. For example, the models A, B, and 
C diffcr subtantially in the mass of their tails despite their 
equal disk-to-halo mass ratio. However, if we use instead the 
ratio of the total amount of mass in the region of the disk to 
the total mass of the galaxy, the mass ratio criterion can be 
partly resurrected. This is shown in the right panel of Fig. 
11, where we plot T eS versus M(2R d )/M. Here M(2R d ) is 
the total mass inside two disk scale lengths, and M is the 
total mass of the galaxy and its halo. In this formulation, 
the mass ratio measures the relative distribution of mass 
within the system, and is a fair indicator of the ability of a 
galaxy model to form tidal tails. However, a detailed com- 
parison with Fig. 10 shows, that £ does a better job than 
M(2R d )/M. For example, the model Tl appears as an out- 
lyer in Fig. 11, failing to fit the monotonic trend of larger 
T e s with increasing M(2R d )/M, but fits within the general 
distribution in Fig. 10. 

Starting from a head on collision, Figure 10 also shows 
that the tidal response T e g becomes larger as the impact 
parameter 7? kep is increased. However, for very wide encoun- 
ters one expects only a small distortion of the disks. As the 
sequence of models C0-C5 demonstrates, there is indeed a 
maximum response for an intermediate impact parameter 
i?ke P of order a few disk scale lengths. 



Also, the orientation of the disks is an important fac- 
tor in determining the strength of thc disk response. In the 
retrograde collision C2r, the spins of the galaxies are just 
reversed compared to C2, yet this already makes the tidal 
tails much weaker, as seen in Figure 12. Similarly, the in- 
clined galaxies of simulation C2i produce tails that are less 
extended than those of run C2. 

We also note that the simulations Tl, Ul, and VI, in 
which our 'disk-dominated' models collide, produce tidal 
tails well in line with the general trend of Figure 10. 

5.4 A model with a bulge 

In Figure 13 we show the time evolution of run Wl, which 
collides two disk + bulge galaxies. These disk galaxies (our 
model W) are descendents of model D, but one third of 
the stellar mass has been put into a centrally concentrated 
bulge. This results in a rotation curve (Figure 3) that is 
practically flat up to the very center of the disk, with a 
shape quite similar to the models of DMH. 

Due to their strong central concentration, the bulges 
survive largely unaffected until their final coalescence. How- 
ever, the disks develop prominent tails that appear to be 
similar in strength to the other models with specinc angular 
momentum corresponding to A = 0.05. A measurement of 
T e ff shows that the strength of the tails is in fact quite sim- 
ilar to the directly comparable simulation Dl. Also, the run 
Wl fits well into the plot of Figure 10, although the inner 
structure of model W is very different from that of the other 
models. 

This suggests that bulges are not effective in preventing 
tail formation, at least as long as they primarily affect the 
inner rotation curve. 



6 DISCUSSION 

In this study we constructed N-body models of disk galax- 
ies with structural properties directly motivated by current 
theories of hierarchical galaxy formation. In particular, thc 
mass of the dark haloes in these models is much larger than 
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Figure 13. Time evolution of run Wl, which collidcs our disk+bulgc model with itself. The top series of pancls displays thc disk particlcs 
projected onto the orbital plane, while the bottom pancls show thc bulge particlcs. Thc lcngth units refer to /i _1 kpc, and the elapsed 
time since the start of the simulation (upper left corners) is given in units of 0.1 Hubble times, or 9.8 X 10 8 7i -1 yr. 
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that of the stellar disks. In the most extreme models we 
consider, the halo-to-disk mass ratio is larger than 40:1. 

Provided the spin parameter is not too small, these 
models can produce long and massive tidal tails, despite 
the massive haloes. Halo-to-disk mass ratio is not a useml 
indicator for tail-making ability. 

Instead, the size of the disk compared to that of the 
halo seems to be the critical factor. In our approach, the 
size of the disk is tied to the spin of the dark halo. A larger 
spin parameter A leads to larger disks. The bulk of the disk 
material is then more loosely bound in the dark matter po- 
tential well and can be more easily induced to form long 
tidal tails. This effect can be quantified in terms of the ratio 
£ of the circular velocity to the escape speed at a radius 
R — 2Rd- We have shown that £ correlates well with the 
tidal response of the disk models. For models with £ > 8 wc 
do not expect signihcant tails, while models with £ < 6.5 
can produce substantial tails. 

When £ is used to characterize the tail-making ability, 
the results of DMH agree with our own. In their sequence 
of four models not only the mass of the halo changes, but 
also its spatial extent. We think the latter effect is criti- 
cal in denning tail-making ability, because the relative size 
of disk and halo affects the value of £ strongly. This also 
agrees with the earlier conclusion of Barnes (1997, private 
communication) . For given inner structure and given mass 
ratio Md/Mh, less extensive halos have larger £ and so make 
weaker tails. 

We focused in this work on just one halo mass. Note 
however, that the shape of the rotation curves does not de- 
pend on our particular choice for «200- Rotation curves with 
other peak velocities may be realized by an appropriate scal- 
ing of V2oo ■ 

According to the work of NFW, the shape of the dark 
matter prohle is insensitive to cosmology. Also, the distri- 
bution of A is universal and independent of the initial power 
spectrum. Furthermore, the average value of the concentra- 
tion c does not vary with halo mass, and it depends only 
weakly on cosmology. For low-density nat universes a smaller 
value, c ~ 5, is probably more appropriate than the value 
c = 15 employed in most of the models in this work. Such 
a smaller concentration is also supported by observational 
data (Navarro 1998). However, our model 'VI', which has 
c = 5, produces tails well in line with the other models in 
Figure 10. Also note that according to Figure 2, smaller c 
gives rise to larger disks, in principle favouring even stronger 
tails. 

This suggests that all 'reasonable' CDM cosmologies 
can produce disk galaxies with A > 0.05 that are roughly 
equally capable of producing tidal tails when they collide 
and merge with similar objects. We conclude that the ob- 
served lengths of tidal tails in interacting galaxies are con- 
sistent with current CDM cosmologies, and that tidal tails 
are not useful to discriminate between different fiavours of 
these scenarios. 

We hope that the N-body representations of disk models 
constructed in this work may be more realistic caricatures 
of real spiral galaxies than those of previous work. In partic- 
ular, the structural properties of our models are motivated 
by hierarchical structure formation and are less ad hoc than 
in previous simulations of this kind. These models should be 
useful for future work on galaxy evolution. 
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